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ABSTRACT 


The  stability  of  the  viscous  flow  between  two  parallel 
horizontal  plates  due  to  a constant  reduced  pressure  gradient 
in  a system  rotating  about  a vertical  axis  is  studied.  The 
critical  value  of  the  Reynolds  number  R,  based  on  the  reduced 
pressure  gradient,  is  a function  of  a dimensionless  rotation 
parameter  T,  the  Taylor  number.  A numerical  solution  of  the 
eigenvalue  problem  shows  that  (i)  the  viscous  instability  mode 
associated  with  plane  Poiseuille  flow  at  T = 0 disappears  at 
a value  of  T of  about  T = 0.4,  and  (ii)  for  T f 0 a new  in- 
stability  mode  appears  as  a result  of  the  Coriolis  effect  on 
the  basic  flow  and  in  the  perturbation  equations.  This  new 
instability  gives  the  critical  value  of  R for  values  of  T as 
small  as  0.06. 


I INTRODUCTION 


The  Navier-Stokes  equations  have  an  exact  solution  for  the 
flow  between  two  horizontal  parallel  planes  due  to  a pressure 
gradient  in  the  presence  of  rotation  about  an  axis  perpendicular 
to  the  planes.  The  rotation  forces  a component  of  flow  in  a 
horizontal  direction  perpendicular  to  the  direction  of  the 
pressure  gradient  as  well  as  a modification  of  the  parallel  com- 
ponent. Wollkind  and  DiPrima^  studied  the  stability  of  this 
flow  for  small  values  of  the  dimensionless  rotation  parameter  T 
[see  Eq.  (1) ] by  means  of  an  expansion  in  T.  In  this  paper  we 
use  direct  numerical  procedures  to  solve  the  stability  problem 
for  moderate  values  of  T. 

The  results  confirm  those  of  Wollkind  and  DiPrima  for  T -►  0 
for  the  instability  of  the  critical  mode  associated  with  the 
stability  of  plane  Poiseuille  for  (T  = 0) . The  calculations 
also  show  the  rather  unexpected  result  (and  one  not  obtainable 
by  the  expansion  procedures  used  by  Wollkind  and  DiPrima)  that 
this  viscous-driven  instability  mode  disappears  at  a value  of  T 
slightly  greater  than  0.35.  Moreover,  they  reveal  the  presence 
of  a second  type  of  instability  when  T ^ 0.  This  new  instability 
gives  the  critical  Reynolds  number  of  the  basic  flow  for  values 
of  T greater  than  about  0.06.  The  critical  disturbance  is  a wave 
propagating  very  slowly  (almost  a standing  wave)  at  an  angle 
(nearly  perpendicular  for  T very  small)  to  the  direction  of 
the  pressure  gradient. 
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The  stability  problem  is  formulated  in  Sec.  II  and  the 
nxamerical  procedure  for  solving  the  eigenvalue  problem  is 
explained  in  Sec.  III.  Results  given  in  Sec.  IV  are  summarized 
in  Sec.  V. 


II  THE  STABILITY  PROBLEM 


The  reader  is  referred  to  Wollkind  and  DiPrima^  for  the 
details  of  the  formulation  of  the  stability  problem.  We  present 
here  only  a summary  sufficient  for  us  to  state  the  stability 
problem. 

Consider  an  xyz  coordinate  system  which  is  rotating  about 
the  z axis  with  constant  angular  velocity  Jl.  We  suppose  that 
there  are  planes  at  z « 0 and  z = d,  and  that  there  is  a constant 
reduced  pressure  gradient  3p/9x  in  the  direction  of  the  positive  x axis. 
Let 


d^  ^ 
8pv  3x  ' 


(1) 


where  v and  P are  the  kinematic  viscosity  and  density,  respectively. 

The  parameter  T,  the  Taylor  nxunber,  is  a measure  of  the  ratio  of 
the  Coriolis  force  to  the  viscous  force.  It  is  the  recipient  of 
the  Ekman  number.  The  pareuneter  U*  is  an  appropriate  velocity 
scale  for  the  case  3p/3x  • constant  (assumed  negative  so  that 
U*  > 0),  and  R is  the  Reynold  number. 

An  exact  solution  of  the  equations  of  motion  satisfying 
the  no  slip  boundary  conditions  at  the  walls  is 
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[-sinh  c sin  C A cosh  c sin  C B slnh  ( cos  c). 
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V_  = T ^ [cosh  C cos  C - 1 - A sinh  C cos  C + B cosh  ? sin  5],  (3) 

A = [sinh  (2T^/^) ]/tcosh  + cos  (2T^/^)]  , 

(4) 

B = (sin  (2T^/^) ]/[cosh  (2t^/^)  + cos  (2T^/^)] 

1/2 

where  c = T z and  z is  now  a dimensionless  variable  scaled  with 
respect  to  d/2  so  that  0 £ z £ 2.  The  dimensionless  velocities 
U_  and  V_  (with  respect  to  U*)  in  the  x and  y directions,  respec- 
tively,  are  symmetric  about  the  center  of  the  channel  (z  = 1) . 

For  T small  the  velocity  reduces  to  a small  perturbation  from 
plane  Poiseuille  flow;  Ug  = z(2-z)  + O(T^)  and  Vg  = T(-z^  + 4z^  -8z)/6 
+ O(T^) . For  T large  and  z bounded  away  from  the  wails,  we  obtain 
the  geostrophic  flow  U_  *'»  0 and  V_  'v*  (3p/3x)/2pnu*  =»  -T  We 
note,  since  3p/3x  < 0 so  that  U*  > 0,  that  the  cross  flow  is  in 
the  negative  y direction  when  the  rotation  is  counterclockwise. 

Next  we  assume  that  all  quantities  have  been  made  dimension- 
less: lengths  with  respect  to  d/2,  velocities  with  respect  to  U*, 
and  time  with  respect  to  d/2U*.  To  study  the  stability  of  the 
basic  flow  (Ug^Vg)  we  superimpose  a small  disturbance,  (i)  sub- 
stitute in  the  equations  of  motion  and  neglect  quadratic  terms, 

(li)  look  for  normal  mode  solutions  of  the  form 

u(x,y,z,t)  - Ug(z)  + u'(z)  exp[i(ax  + Py  - kct) 1 , (5) 

2 2 1/2 

where  k ■ (o  + 6 ) » (iii)  introduce  new  coordinates  (Xj^,yj^) 

with  velocity  perturbation  eunplitudes  u^^  and  Vy  respectively,  by 
rotating  the  (x,y)  coordinate  system  through  an  angle  6 measured 
counterclockwise  from  the  positive  x direction  so  that  Xj^  « 

X cos  6 y sin  6,  * -x  sin  6 4-  y cos  6,  where  cos  6 > a/k 
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and  sin  0 = B/k,  (iv)  multiply  the  x-momentum  equation  by  (-o) , 
the  y-momentum  equation  by  0 and  add  to  obtain  a new  equation, 

(v)  multiply  the  x-momentum  equation  by  (-6),  the  y-momentum 
equation  by  a and  add  to  obtain  a second  equation,  and  (vi)  eliminate 
the  pressure  by  using  the  z-momentum  equation,  and  use  the  continuity 
equation  to  introduce  a function  41^  such  that  Uj^  = -d(|)j^/dz  and 
w'  = ik((»j^. 

The  result  of  these  calculations  is  the  following  sixth  order 
system  of  ordinary  differential  equations: 

(D^  - - ikR[(Uj^-c)  (D^-k^)(()j^  - (D^Uj^)(j)j^]  - 2TDVj^  = 0,  (6) 


(D^  - k^)Vj^  - ikR[  (Uj^-c)Vj^  + + 2TD(j)j^  = 0 . 


The  velocities  and  are  given  by 

Uj^(2,T,0)  = Ug(z,T)  cos  0 + Vg(z,T)sin  0 , 

Vj^(2,T,0)  a -Ug(z,T)sin  0 + Vg(z,T)cos  0 , 

and  in  the  (Xj^,yj^,z)  coordinate  system  the  perturbation  is  pro- 
portional to  exp  lik(Xj^-ct)  ] . 

The  boundary  conditions  are 

♦ 1 a D(>J^  a a 0 at  Z a Q , 

and  for  disturbances  for  which  is  symmetric  and  v^  is  anti- 
symmetric about  z a 1, 

a a Vj^  a 0 at  2 a 1. 
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For  disturbances  for  which  (|)j^  is  anti-synunetric  and  is  syiranetric 

2 

about  z = 1,  the  boundary  conditions  would  be  = D = 0 

at  z = 1.  It  is  known  for  T = 0 that  the  critical  disturbance  is 
one  for  which  is  synunetric  about  z = 1,  and  only  such  dis- 
turbances will  be  considered  here. 

Equations  (6),  (7),  (10),  and  (11)  define  an  eigenvalue 

problem  of  the  form  F(T,R,9,k,c)  = 0.  The  Taylor  number  T,  the 

Reynolds  nujnber  R,  and  the  angle  of  propagation  9 are  real.  For 

temporally  growing  or  decaying  disturbances,  the  wavenumber  k is 

real  and  the  parameter  c (the  complex  cunplif ication  rate)  is  complex, 

c = c^  + ic^.  For  spatially  growing  or  decaying  disturbances  k is 

complex  with  kc  real.  In  this  paper  we  consider  the  case  of 

temporal  instability.  Thus,  k is  real  and  the  flow  is  unstable 

if  there  exists  an  eigenvalue  c such  that  c^  > 0.  For  a given 

value  of  T the  condition  that  the  eigenvalue  with  largest  imaginary 

part  have  c^  = 0 determines  a neutral  surface  in  R,9,k  space  which 

separates  stable  from  unstable  states.  The  critical  value  of  R 

is  the  smallest  value  of  R that  corresponds  to  a point  on  the 

neutral  surface:  R (T)  = min  R(T,k, 9 ,c^,c . * 0)  for  k > 0 

c r 1 

-Tr/2  < 9 i ■it/2.  (Note  that  0 can  be  restricted  to  this  interval 
since  c^  can  be  positive  or  negative.)  The  corresponding  values 
of  k,0,c^  determine  the  wavenumber,  direction  of  propagation, 
and  the  wave  velocity  of  the  critical  disturbance. 

Before  turning  to  the  numerical  procedure  for  solving  the 
eigenvalue  problem,  it  is  useful  to  derive  an  "energy"  integral 
associated  with  Eqs.  (6)  and  (7).  Multiplying  Eqs.  (6)  and  (7) 
by  the  complex  conjugates  of  and  respectively.  Integrating 


6. 


by  parts,  adding,  and  taking  the  real  part,  we  obtain 
kRc. (E.  + E„)  = -(D.  + D„)  + R(M  + H)  . 

1 (J>  V (})  V 


(12) 


Here 


1 

0 


1 

Dv  = y*  + k^lv^l^)dz  , (13) 

0 


1 1 

M = k y (DUj^)  (().j^j.D(J)3^^-(J.^.D<^^^)dz  , H = ky*  . 

0 0 


The  quantity  kRc^ (E^  + E^)  can  be  interpreted  as  the  rate  of  change 
of  disturbance  kinetic  energy,  D.  + D is  the  rate  of  disturbance 

9 V 

viscous  dissipation,  RM  is  the  rate  of  transfer  of  kinetic  energy 
from  the  component  of  the  basic  velocity  to  the  disturbance, 
and  RH  is  the  rate  of  transfer  of  kinetic  energy  from  the  Vj^ 
component  of  the  basic  velocity  to  the  disturbance. 

Since  E,,  E , D . , D„  are  positive  definite,  it  follows  that 

(p  V <p  V 

M and/or  H must  be  positive  if  disturbances  are  to  grow  (c^^  > 0)  . 

If  M is  positive,  disturbance  energy  is  gained  from  the  Uj^  component 
of  the  basic  velocity  which  is  parallel  to  the  line  of  propagation 
of  the  disturbance,  while  if  H is  positive,  disturbance  energy  is 
gained  from  the  Vj^  component  of  the  basic  velocity  which  is  per- 
pendicular to  the  line  of  propagation  of  the  disturbance. 


It  is 
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interesting  that  T does  not  appear  explicitly  in  Eq.  (12) ; it 
enters  only  through  its  effect  on  the  basic  velocity  field  and 
the  solution  of  the  perturbation  equations. 


Ill  NUMERICAL  PROCEDURE 

It  is  well  known  that  numerical  solutions  of  Orr-Sommerfeld 

like  problems  are  difficult  to  obtain,  because  of  the  presence  of 

solutions  that  grow  very  rapidly  in  z when  the  Reynolds  number 

is  large.  The  difficulties  can  largely  be  overcome  by  using 

2 3 

orthogonalization  techniques  (cf.,  Conte  , Davey  , and  Scott  and 

4 

Watts  ).  For  the  problem  under  consideration,  we  followed  Davey 's 
method,  which  uses  parallel  shooting  and  orthonormalizations  to 
solve  the  system  (6)  to  (11)  for  prescribed  values  of  T,  R,  6,  k 
and  an  initial  guess  for  the  eigenvalue  c.  Successive  iterates 
for  c were  obtained  using  Muller’s  method^  to  find  a root  of  the 
transcendental  equation  F(T,R,6,k,c)  = 0. 

The  differential  equations  (6)  and  (7)  were  integrated  using 
an  implicit  Adams  method  with  a step  size  of  1/250  and  four 
orthonormalizations  for  Reynolds  numbers  of  approximately  6000. 

In  order  to  estimate  the  accuracy  of  the  numerical  solutions, 
we  (i)  calculated  several  eigenvalues  for  plane  Poiseuille  flow 
and  compared  results  with  Davey^,  (ii)  calculated  solutions  of  a 
model  sixth  order  problem  having  a known  exact  solution,  and 
(iii)  calculated  residuals  of  the  differential  equations  (6)  and  (7) 
using  the  computed  values  of  Vj^,  and  c.  In  general,  we 
estimate  that  the  eigenvalues  and  eigenfunctions  are  accurate 
to  six  or  seven  significant  digits. 
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Our  primary  goal  was  to  find  points  on  the  neutral  surface. 
Thus,  once  a satisfactory  eigenvalue  c was  found  by  Davey's 
procedure,  we  changed  R and/or  0 and  used  linear  extrapolation 
to  search  for  a point  on  the  neutral  surface.  An  attempt  was  made 
to  change  R and/or  0 so  that  the  search  direction  was  approximately 
orthogonal  to  the  neutral  surface. 

Once  a point  on  the  neutral  surface  was  found,  the  corres- 
ponding eigenvector  (4)j^,Vj^)  was  calculated  using  a technique 

2 4 

similar  to  those  used  by  Conte  and  Scott  and  Watts  . The  entire 
process  worked  satisfactorily  provided  a sufficiently  close  initial 
guess  for  a point  on  the  neutral  surface  was  known.  This  was 
generally  easy  to  accomplish  by  continuation  from  results  already 
obtained.  Typically,  four  to  six  Muller  iterations  were  required 
to  find  an  eigenvalue  and  three  or  four  linear  extrapolations 
were  required  to  find  a point  on  the  neutral  surface;  thus,  the 
differential  equations  (6)  and  (7)  had  to  be  solved  from  12  to  24 
times  for  each  point  on  the  neutral  surface. 

Finally,  minimum  points  on  the  neutral  surface  were  obtained 
by  adjusting  k and/or  0 using  a pattern  search  procedure  (cf., 
Jacoby,  Kowalik,  and  Pizzo^,  Chapter  4)  to  locate  the  minimum 
value  of  R for  a fixed  value  of  T.  This  procedure  also  worked 
quite  well  provided  that  relatively  small  search  steps  were  used. 


9. 


TV  RESULTS 

7 

For  T ~ 0 it  is  known  (Squire's  Theorem  ) that  the  critical 
Reynolds  number  for  the  eigenvalue  problem  (6) -(11)  occurs  for 

g 

0=0.  The  values  given  by  Reynolds  and  Potter  are 

0=0,  k = 1.0207  , R = 5772.12  , c = 0.26402.  (14) 

c c c ' rc 

Wollkind  and  DiPrima^  found  for  small  T that  this  critical  point 
was  perturbed  as  follows: 

0 = -2.044T  + O(T^)  , 

c 

k = 1.0207  + O(T^)  , 

(15) 

R^  = 5772.12(1  - 0.3877T^  + 0(T^)]  , 

c^^  = 0.26403  + O(T^)  . 

They  noted  that  these  results  could  only  be  expected  to  be  valid 
for  very  small  T and  suggested  that  they  are  probably  not  valid 
at  T = 0.5. 

We  used  the  method  described  in  Sec.  Ill  to  obtain  k , 0 , R , 

c c c 

and  c for  T = 0.05,  0.1,  0.2,  0.3,  0.4,  0.5,  0.75,  and  1.0.  The 

]CC 

results  are  given  in  Table  1.  Observe  that  for  T = 0.05,  0.1,  0.2, 
and  0.3  there  are  two  relative  minima:  one  near  6=0  and  k = 1.02 
and  one  near  0 = 1.3  to  1.5  and  k = 1.37.  The  first  of  these,  which 
we  denote  by  (I) , is  clearly  a small  perturbation  due  to  rotation  of 
the  critical  parameters  for  the  instability  of  plane  Poiseuille 
flow  (T  » 0) . The  perturbation  formulas  (15)  give  results 
in  reasonable  agreement  with  the  numerical  calculations 
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of  the  critical  values  for  the  Type  I disturbance.  The  comparison 
is  shown  in  Table  2.  The  second  minimum,  which  we  denote  by 
(II) , is  a new  instability  introduced  by  the  rotation;  it  will  be 
discussed  later  in  this  section.  We  note  that  two  relative 

9 

minima  on  a neutral  curve  were  found  by  Gill  and  Davey  in  their 
analysis  of  the  instability  of  a buoyancy  driven  system,  and  by 
Lilly^^  in  his  analysis  of  the  instability  of  the  Ekman  boundary 
layer.  Each  of  these  eigenvalue  problems  involves  a sixth  order 
system  of  differential  equations  that  are  similar  to  Eqs.  (6,'  and 
(7);  however,  the  domain  for  each  problem  is  (0,®)  in  contrast  to 
the  finite  domain  of  the  present  problem. 

The  two  relative  minima  of  the  neutral  surface  are  shown  in 
Fig.  1 where  the  variation  of  R on  the  neutral  surface  is  given 
as  a function  of  0 for  k = 1.02  for  the  cases  T = 0,  0.05,  and  0.2. 

Similar  curves  are  given  in  Fig.  2 for  T = 0.3  and  0.5.  The 
minimum  point  corresponding  to  the  type  I instability  disappears 
at  a value  of  T between  0.3  and  0.5  as  can  be  seen  from  Fig.  2 
(it  was  not  possible  to  locate  a minimum  at  T = 0.4).  The  width 
of  the  valley  associated  with  the  type  II  instability  is  extremely 
narrow  for  very  small  values  of  T,  but  increases  as  T increases. 

The  curves  in  Figs.  1 and  2 are  incomplete.  While  we  were  able 
to  take  our  calculations  up  to  Reynolds  numbers  as  large  as  500,000 
and  obtain  (0,R)  points  not  shown  in  Figs.  1 and  2,  we  were  not 
able  to  determine  if  the  curves  have  a finite  maximum  or  go  to 
infinity  (as  is  the  case  for  T = 0,  0 -*•  + 11/2)  . 

The  corresponding  variation  of  with  0 for  k “ 1.02  and  T = 0, 
0.5,  and  0.2  is  shown  in  Fig.  3.  We  note  that  since  R is  a periodic 
function  of  0 with  period  it  it  follows  that  c^(0=it/2)  “ - c^(0=-it/2). 
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It  is  clear  from  Table  1 and  Fig.  1 that  the  absolute  minimum 
on  the  neutral  surface  as  T ->  0 corresponds  to  a type  I instability. 
However,  the  type  II  instability  rapidly  becomes  dominant  and  gives 
the  absolute  minimum  at  a value  of  T as  small  as  0.06  (approximately). 
The  variation  of  R , c , 6 , and  k with  T is  shown  in  Figs.  4,  5, 

6,  and  7.  The  variation  of  with  T is,  of  course,  continuous; 

however,  each  of  the  other  parameters  change  discontinuously  when 
the  type  I instability  is  replaced  by  the  type  II  instability. 

The  eigenfunctions  corresponding  to  the  parameter  values  for 
each  of  the  two  relative  minima  for  T = 0.05  are  shown  in  Figs.  8 
and  9.  The  normalization  for  these  eigenfunctions  as  well  as  the 
case  shown  in  Fig.  10  is  <)>j^(l)  = 1.  The  (|)j^  component  of  the  eigen- 
function for  the  type  I instability  for  T = 0.05  (Fig.  8a  and  8b) 
is  very  similar  to  the  corresponding  eigenfunction  for  the  insta- 
bility of  plane  Poiseuille  flow  (T  = 0);  for  example,  see  Fig.  3a 
of  Reynolds  and  Potter  , and  for  a larger  value  of  R see  Fig.  la 

of  Lee  and  Reynolds^^.  Note  the  relative  simplicity  and  the 

different  order  of  magnitudes  of  the  eigenfunction  for  the  type  II 
instability  (especially  the  Vj^  component)  compared  to  the  eigen- 
function for  the  type  I instability.  The  eigenfunction  associated 
with  the  critical  parameter  values  on  the  neutral  surface  for 
T = 1.0  is  shown  in  Fig.  10.  The  simplicity  of  the  type  II 
instability  eigenfunction  suggests  the  instability  might  be 
calculated  by  procedures  such  as  the  Galerkin  method;  however, 
we  have  not  pursued  this  idea. 

We  turn  now  to  the  type  II  instability  which  determines  the 
critical  Reynolds  number  for  T > 0.06  (approximately).  Notice 
from  Table  1 that  the  value  of  is  extremely  small  and  positive; 
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hence  the  disturbance  is  almost  a standing  wave  moving  very  slowly 

in  the  direction  given  by  6^.  For  small  T we  have  0^  = 7t/2;  hence 

the  direction  of  propagation  is  nearly  perpendicular  (positive  y 

direction)  to  the  direction  of  the  basic  flow  when  T = 0 (the 

positive  X direction) . Moreover,  we  note  that  the  positive  y 

direction  is  opposite  to  the  direction  of  the  geostrophic  flow 

(negative  y direction)  that  occurs  when  T is  large.  With  increasing 

T the  angle  of  the  direction  of  propagation  decreases  from  it/2, 

and  at  T = 1 we  have  0 =0.86  rad  = 50“ . 

c 

There  are  two  ways  in  which  the  Coriolis  effect  manifests 
itself  in  the  stability  analysis.  Firstly,  it  modifies  the  basic 
velocities  and  given  by  Eqs.  (8),  (9),  (2),  and  (3)  that 
appear  in  the  distrubance  equations  (6)  and  (7) ; secondly,  it 
couples  Eqs.  (6)  and  (7)  through  the  terms  -2TDVj^  in  the  equation 
for  and  the  term  2TD())j^  in  the  equation  for  v^^. 

It  is  well  )<nown  that  the  eigenvalue  problem  for  the  Orr- 

Sommerfeld  equations  [Eq.  (6)  with  T = 0)  is  very  sensitive  to 

the  basic  velocity  profile  and  that  inflection  points  play  a 

special  role.  For  T = 0.2  the  velocity  profiles  Uj^(z,0)  are 

shown  in  Fig.  11  for  several  values  of  0 in  the  neighborhood  of 

0 = 1.404.  It  is  clear  from  Fig.  11  that  U,  (z,0)  has  an  inflection 

c i 

point  for  0 = 1.404  and  for  nearby  values  of  0.  Indeed,  a simple 
analysis  shows  that  Uj^(z,0)  first  exhibits  an  inflection  point 
at  z = 1 when  0 = tan”^(cot/T  coth/f ) , and  that  the  inflection 
point  decreases  from  z = 1 to  z = 0 where  it  disappears  when  0 has 
increased  to  it/2. 
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The  interval  in  0 for  which  Uj^(2,G)  has  an  inflection  point 

corresponds  almost  exactly  to  the  valley  near  0 = ti/2  in  the  R 

vs.  0 curve  shown  in  Fig.  1.  For  other  small  values  of  T the  situation 

is  the  same;  the  velocity  profile  Uj^(2,0)  has  an  inflection  point 

for  each  value  of  0 in  a small  interval  about  the  corresponding 

value  of  0^  for  the  type  II  instability. 

This  suggests  that  the  type  II  instability  is  an  inviscid 

instability  associated  with  the  inflection  point  that  occurs  in 

the  component  of  the  basic  velocity  profile  as  a result  of  the 

Coriolis  force.  [For  T = 0 we  have  = z(2-z)cos  0 so  the 

inflection  point  is  a Coriolis  effect.]  In  Table  1 we  have  given 

the  value  of  z*  at  which  the  inflection  point  occurs  for  0 and 

c 

the  value  of  U,  (z*,0  ).  It  is  known  if  viscosity  is  neglected  in 
X c 

the  Orr-Somroerf eld  equation  (R  ->  <*>)  that  for  neutral  disturbances 

Ul  - c^  must  vanish  at  least  once  for  0 < z < 1.  Such  points  are 

critical  points  of  the  inviscid  disturbance  equation.  For  velocity 

profiles  with  a point  of  inflection  z*  and  only  one  critical  point, 

we  have  c^  equal  to  the  velocity  at  the  inflection  point.  However, 

for  velocity  profiles  with  an  inflection  point  and  two  critical 

points,  which  is  typical  of  several  of  the  velocity  profiles  shown 

12 

in  Fig.  11,  Gregory,  Stuart,  and  V'Jalker  (see  Section  12)  have 
shown  that  c^  will  be  slightly  less  than  the  velocity  at  the  in- 
flection point.  It  is  clear  from  Table  1 that  Uj^(z*,0^)  and  c^ 
are  of  the  same  (small)  size  with  = 2U,  (z*,0^).  It  is  significant 
that  c and  U, (z*,0  ) have  comparable  values.  We  have  no  explanation 

IT  X O 

as  to  why  c » 2U, (z*,0  ) rather  than  c ^ U, (z*,0  );  except  to  note 
that  viscous  and  Coriolis  effects  are  not  included  in  this  discussion. 


1.4  . 

The  importance  of  the  modification  of  the  component  of  the 
basic  velocity  due  to  rotation  can  be  shown  as  follows.  Firstly, 
consider  Eq.  (6)  with  T = 0,  the  Or r-Sommorf eld  equation.  Then 
= z(2-z)cos  0 and  at  k = 1.02,  0 = 1.4  3 the  critical  value  of 
R is  5772.12  seed. 43)=  41,132.  Now  for  T = 0.2  we  neglect  the 
Coriolis  term  (-2TDVj^)  in  Eq.  (6),  but  use  the  velocity  profile 
that  is  correct  for  this  value  of  T.  Then  we  find  that  for  k = 1.02 
the  minimum  point  on  the  neutral  curve  occurs  for  0 = 1.43  and  has 
the  value  R = 6688.  This  is  a significant  change  from  R ^ 41,000. 
Moreover,  the  Orr-Sommerfeld  equation  with  the  Uj^  velocity  profile 
corresponding  to  T = 0.2  gives  the  qualitative  features  of  the 
full  sixth  order  system  in  the  neighborhood  of  the  type  II  minimum 
point  as  can  be  seen  from  Fig.  12.  On  the  other  hand,  modification 
of  the  basic  velocity  by  rotation  is  not  the  entire  answer. 

For  the  full  sixth  order  system,  Eqs.  (6)  and  (7) , the  minimum 
point  on  the  neutral  surface  for  T = 0.2,  k = 1.02  is  R = 1953 
at  0 = 1.40. 

Th:*  energy  integral  relation  (12)  is  also  useful.  In  Table  3 
the  values  of  the  different  integrals  have  been  calculated  for  the 
minimum  points  on  the  neutral  curves  corresponding  to  type  I and 
type  II  disturbances.  In  these  calculations,  the  eigenfunctions 
have  been  further  normalized  by  the  requirement  = 1.  Also, 
the  integrals  were  evaluated  using  the  trapezoidal  rule  with  an 
interval  size  of  1/50  for  most  of  the  calculations;  hence  the 
evaluation  of  the  integrals  is  not  as  accurate  as  the  eigenfunctions. 
Recall  that  RM  and  RH  represent  the  rate  of  transfer  of  kinetic 
energy  from  the  Uj^  and  the  Vj^  components  of  the  basic  velocity. 
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respectively,  to  the  disturbance.  Consider  the  type  l instability. 
For  T = 0.05,  when  the  type  I disturbance  is  the  critical  dis- 
turbance, the  primary  contribution  to  the  disturbance  energy 
is  from  the  components  as  would  be  expected.  However,  with 
increasing  T the  situation  changes  rapidly  and  et  T = 0.3  the 
primary  contribution  (by  an  order  of  magnitude)  comes  from  the 
perpendicular  component  of  the  basic  velocity. 

For  the  type  II  instability,  essentially  all  of  the  dis- 
turbance energy  is  contributed  by  the  interaction  of  with  <))j^ 
and  v^ . It  is  clear  from  Table  3 that  the  energy  balance  for 
a type  li  neutral  disturbance  is  0 = -D^  + RH.  Thus,  the  type  II 
instability  appears  to  be  driven  by  the  Coriolis  force  through 

(i)  modification  of  the  basic  component  of  velocity,  and 

(ii)  the  transfer  of  energy  through  the  interaction  of  the 
Reynolds  stress  and  the  gradient  of  the  component  of  velocity. 

Finally,  we  consider  the  direct  coupling  effect  of  the 
Coriolis  force  in  Eqs.  (6)  and  (7).  When  T = 0,  Eq.  (6)  with 
the  boundary  conditions  on  <|>j^  from  Eqs.  (10)  and  (11)  can  be 
used  to  determine  the  eigenvalue  relation.  Then  the  non- 
homogeneous  differential  equation  (7)  can  be  solved  for  Vj^ 
subject  to  the  boundary  conditions  (10)  and  (11),  provided 
that  the  corresponding  homogeneous  problem  does  not  have  a 
solution.  It  can  readily  be  demonstrated  as  first  shown  by 
Squire^  that  the  latter  problem  only  has  solutions  corresponding 
to  decaying  disturbances  (c^  < 0) ; hence,  no  difficulty  will  be 
encountered  in  calculating  neutral  or  amplified  disturbances. 
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(Nor  is  it  likely  that  an  eigenvalue  < 0 of  the  boundary  value 
problem  for  will  be  an  eigenvalue  of  the  homogeneous  boundary 
value  problem  for  Indeed,  this  calculation  forms  the  first 

term  in  a perturbation  calculation  in  powers  of  T (see  Wollkind 
and  DiPrima^) . 

However,  there  is  another  solution  when  T = 0;  namely,  the 
solution  = 0 with  the  eigenvalue  relation  determined  by 

(D^  - k^)v^  - ikR(U^  - c)v^  = 0 , v^(0)  = v^(l)  = 0.  (16) 

As  noted  above,  it  is  not  difficult  to  show  that  the  eigenvalues  (c) 
of  this  problem  give  decaying  disturbances.  Indeed,  if  we  write 
Uj^(z,e)  = Uj^q(z)cos  0,  c = CqCOs  0,  and  R » Rq/cos  0,  where 
Uio(z)  * z(2-z),  and  then  multiply  the  equation  by  the  complex 
conjugate  of  Vj^  and  integrate,  we  find  0 < < 1 and 


c . » — ^ ^ 

Oi  kRo 


y*(|DVj^|^  + k^|Vj^|^)dz 


/iv.i 


- 1 , 2 ^ . 2, 

- -kRr<’  + k ) . 


'dz 


(17) 


Equation  (16)  is  for  an  antisymmetric  solution;  for  a symmetric 

2 2 
solution  the  it  in  Eq.  (17)  is  replaced  by  n /4.  Thus,  the 

symmetric  mode  with  » 0 may  be  less  stable  than  the  anti- 
symmetric mode  with  ■ 0. 

We  have  not  tried  to  study  how  the  eigenvalues  of  the  two 

problems  move  for  T 0.  This  is  an  interesting  mathematical 
problem  whose  analysis  may  clarify  the  origin  of  the  type  II 
instability. 


■s- 


17. 


p'  t V SUMMARY 

The  present  numerical  calculations  show  (i)  the  type  1 
instability  associated  with  the  viscous  instability  of  plane 
Poiseuille  flow  gives  the  critical  Reynolds  number  for  T 0, 

(ii)  the  type  I instability  disappears  at  a Taylor  number  of 
about  0.4,  (iii)  there  is  a type  II  instability  due  to  Coriolis 
effects,  (iv)  the  type  II  instability  gives  the  critical  Reynolds 
number  for  Taylor  numbers  slightly  greater  than  0.06  (approximately), 
and  (v)  the  type  II  instability  represents  a wave  propagating  very 
slowly  (almost  a standing  wave)  along  a direction  at  an  angle 
measured  coxinterclockwise  from  the  direction  of  the  applied 
pressure  gradient  that  decreases  from  nearly  90®  for  T ^ 0.06 
to  about  50®  at  T = 1. 

It  is  interesting  that  an  almost  standing  wave  solution  also 
occurs  at  one  of  the  two  relative  minima  on  the  neutral  surface 
in  the  stability  of  the  E)anan  boundary  layer^®  (T  -*■  *)  . However, 
for  that  problem  the  stationary  wave  solution  does  not  give  the 
absolute  minimum  of  the  neutral  surface.  Rather  stationary  waves 
become  unstable  at  a Reynolds  number  of  about  twice  the  critical 
Reynolds  number  for  the  growth  of  a viscous-Coriolis  disturbance 
which  propagates  at  a small  angle  from  the  direction  of  the 
geostrophic  flow.  It  is  tempting  to  conjecture  that  for  in- 
creasing T the  viscous  instability  mode  of  plane  Poiseuille  flow 
is  replaced  by  an  inflection  point  - Coriolis  instability  at 
T • 0.06  which  persists  as  an  almost  standing  wave  solution  for 
all  T,  but  is  replaced  (in  determining  R ) by  a viscous-Coriolis 
mode  for  large  T. 

k-  w 
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The  destabilizing  effect  of  even  a small  eunount  of  rotation 

(small  T)  is  clearly  shown  in  Fig.  4;  for  example,  the  values  of 

at  T = 0,  0.1,  and  0.2  are  5772,  3575,  and  1806,  respectively, 

with  R = 468  at  T = 1.  It  also  appears  from  Fig.  4 that  the 
c 

curve  of  R^  vs.  T has  nearly  reached  a minimum  at  T = 1,  and 

beyond  this  minimum  R^  will  start  to  increase  with  increasing  T. 

3/2 

This  would  be  consistent  with  the  asymptotic  value  R 55T 

as  T given  by  Wollkind  and  DiPrima^  which  is  based  on 

Lilly' s^^  analysis  of  the  stadiility  of  the  Ekman  boundary  layer. 

Finally,  we  note  that  the  values  of  R^  obtained  here  for 

0 ^ T i 1 lie  above  (as  they  should)  the  energy  limit  for  the 

13 

growth  of  disturbances  calculated  by  Jankowski  and  Squire 
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4 Table  1.  Minimum  points  on  the  neutral  surface 


T 

k 

e 

R 

z* 

Uj^(z*,e 

Type  I Disturbance 

0.05 

1.021 

-0.1026 

5767 

0.2635 

0.1 

1.022 

-0.2061 

5750 

0.2620 

0.2 

1.027 

-0.4234 

5674 

0.2553 

0.3 

1.039 

-0.6786 

5519 

0.2413 

0.35 

1.052 

-0.8577 

5387 

0.2263 

Type  II  Disturbance 

0.04 

1.367 

1.537 

8911 

0.002208 

0.607 

0.00127 

0.05 

1.370 

1.529 

7132 

0.002634 

0.596 

0.00127 

0.1 

1.368 

1.487 

3575 

0.005441 

0.600 

0.00277 

0.2 

1.367 

1.404 

1806 

0.01100 

0.601 

0.00567 

0.3 

1.369 

1.324 

1224 

0.01594 

0.599 

0.00804 

0.4 

1.369 

1.246 

939.5 

0.02081 

0.599 

0.0105 

0.5 

1.3697 

1.1713 

772.9 

0.02533 

0.599 

0.0128 

0.75 

1.3717 

1.0025 

462.6 

0.03517 

0.600 

0.0180 

1.0 

1.3744 

0.85958 

468.5 

0.04279 

0.600 

0.0221 
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Table  2.  Comparison  of  perturbation  calculations  of 
Wollkind  and  DiPrima  (W  & D)  and  present 
numerical  calculations  (F  & D)  for  Type  I 
disturbances. 


Calculation 

T 

®c 

R 

c 

°rc 

W & D 

0.05 

1.0207 

-0.1022 

5767 

0.26403 

F & D 

1.021 

-0.1026 

5767 

0.2635 

W & D 

0.1 

1.0207 

-0.2044 

5750 

0.26403 

F & D 

1.022 

-0.2061 

5750 

0.2620 

H & D 

0.2 

1.0207 

-0.4088 

5683 

0.26403 

F & D 

1.027 

-0.4234 

5674 

0.2553 

W & D 

0.3 

1.0207 

-0.6132 

5571 

0.26403 

F fc  D 

1.039 

-0.6786 

5519 

0.2413 

W & D 

0.35 

1.0207 

-0.7154 

5498 

0.26403 

F « D 

1.052 

-0.8577 

5387 

0.2263 

Table  3.  Energy  integrals  for  critical  values  of  the  parcuneters 


for  Type  I and  Type  II  disturbances  with  = 1.  The 


values  of  R 

c 

' ®c'  ^^c' 

and  c for 

rc 

a given  value  of  T 

are  given  in 

Table  1. 

Numers  in 

parentheses 

indicate 

the  power  of 

ten. 

T 

M 

H 

D.  + D 

(J>  V 

R(M+H) 

Type  I Disturbance 

0.05 

53.95 

2312. 

188.9 

0.3990 

0.03280 

2501. 

2490. 

0.1 

13.17 

560.6 

187.4 

0.09701 

0.03266 

748.1 

745.6 

0.2 

2.899 

119.8 

181.6 

0.02096 

0.03206 

301.3 

300.8 

0.3 

0.9513 

36.98 

169.8 

0.006628 

0.03082 

206.8 

206.7 

Type  II 

Disturbance 

0.04 

0.3475(-4) 

0.3577 (-3) 

11.87 

0.6334(-8) 

0.1331(-2) 

11.87 

11.87 

0.05 

0.5482(-4) 

0.5625(-3) 

11.88 

0.1265(-7) 

0.1665(-2) 

11.88 

11.88 

0.1 

0.2177 (-3) 

0.2239(-2) 

11.87 

0.9926(-7) 

0.3320(-2) 

11.87 

11.87 

0.2 

0.8664(-3) 

0.8926(-2) 

11.86 

0.7778(-7) 

0.6573(-2) 

11.87 

11.87 

0.3 

0.1956(-2) 

0.2012(-1) 

11.87 

0.2607  (-5) 

0.9708(-2) 

11.89 

11.89 

0.5 

0.5420(-2) 

0.5579(-l) 

11.87 

0.1148(-4) 

0.1542(-1) 

11.93 

11.93 

0.75 

0.1213(-1) 

0.1250 

11.88 

0.3546(-4) 

0.2130(-1) 

12.00 

12.00 

1.0 

0.2139(-1) 

0.2208 

11.89 

0.7567(-4) 

0.2577(-l) 

12.11 

12.11 
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Legends  for  the  Figures 

Variation  of  R on  the  neutral  surface  with  0 for  k = 1.02 
and  for  T = 0#  0.05,  and  0.2. 

Variation  of  R on  the  neutral  surface  with  0 for  k = 1.02 
and  for  T = 0.3  and  0.5. 

Variation  of  c^  on  the  neutral  surface  with  0 for  k = 1.02 
and  for  T <=  0,  0.05,  and  0.2. 

Variation  of  R with  T. 
c 

Variation  of  c with  T. 
rc 

Variation  of  0_  with  T. 
c 

Variation  of  k_  with  T. 
c 

Type  I eigenfunction  for  T = 0.05,  0 =■  -0.1026,  k = 1.021, 

R - 5767.  (a)  , (b)  Vj^. 

Type  II  eigenfunction  for  T = 0.05,  0 « 1.529,  k = 1.370, 

R - 7132.  (a)  , (b)  v^. 

Type  II  eigenfunction  for  T > 1.0,  0 « 0.85958,  k » 1.374, 

R - 468.45.  (a)  4^,  (b)  v^. 

U^(z,6)  for  T >■  0.2  and  several  values  of  0. 

Critical  values  of  R as  a function  of  0 for  T * 0.2,  k - 1.02 
Orr-Sonunerfeld  equation;  sixth  order  problem. 
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Poiseuille  flow  at  T = 0 disappears  at  a value  of  T of  about 

T - 0.4,  and  (ii)  for  T 7^  0 a new  instability  mode  appears  as  a 
result  of  the  Coriolis  effect  on  the  basic  flow  and  in  the 
perturbation  equations.  This  new  instability  gives  the  critical 
value  of  R for  values  of  T as  small  as  0.06. 
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